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We calculate the upper critical induction B C 2(8,</>) of an orthorhombic ferromagnetic supercon- 
ductor with parallel-spin p-wave pairing along (=3 and an ellipsoidal Fermi surface with effective 
masses m\,m,2, and m.3, where B ■ e.3 = cos 8. At fixed 8, an ordinary B C 2{<f>) is obtained, but at 
fixed cj>, a novel B C 2(8) is predicted. For fixed 3 < 7713 /(mi cos 2 (f> + m2 sin 2 <f>), an anomalous peak in 
B C 2 occurs at 0° < 8* < 90° , providing a sensitive bulk test of the order parameter orbital symmetry 
in both phases of clean URhGe samples. The method is generalizable to any order parameter. 



Recent discoveries of materials with coexistent super- 
conductivity and ferromagnetism spawned renewed in- 
terest in parallel-spin triplet superconductivity, the sim- 
plest cases having p-wave orbital symmetry [ll4ll | . In a 



non-cubic crystal structure, there can be a variety of dif- 
ferent p-wave states present, which can be characterized 
by measurements of the temperature T dependence of 
the upper critical induction J5 c2 (l2l - fl8j |. Ambient pres- 
sure measurements of the ferromagnetic superconductors 
UCoGe 0, and URhGe 043 showed that the super- 
conductivity exists completely within the ferromagnetic 
T range and that the same electrons are responsible for 
the superconductivity and the ferromagnetism^. 

In URhGe, B C 2 measurements in the low-field regime 
on a sample with room temperature resistivity ratio 
(RRR) = 21 were fit to the Scharnberg-Klemm theory 
of the crystal-aligned p-wave polar state of completely 
broken symmetry (CBS), using the three resistively- 
measured slopes of B C 2 just below the ferromagnetic de- 
magnetization jumps at T c as the only fitting parameters 
The measured B C 2. a {T) fit the predicted behavior 
for the polar state, and both B C 2,b(T) and B C 2 jC (T) fit 
the qualitatively different predicted curve for the CBS 
state, with a constant ratio B C 2^{T) / B C 2 tC {T) due ap- 
parently to T-independent "effective mass" anisotropy. 
In addition, even in the low-field superconducting regime, 
in all three crystal directions, H C 2 strongly violated the 
Pauli limit H p ~ 1.85T C (K) for a singlet-spin s-wave 
superconductor :4j, further suggesting a parallel-spin pair 
state coexisting with the ferromagnetic state. Together, 
the data provided strong evidence that the superconduct- 
ing order parameter is likely to have the simplest parallel- 
spin p-wave orbital form consistent with ferromagnetism 
in the be plane of an orthorhombic crystal, dk a , where 
the pair-spin d = (b + ic)/-\/2, with the p-wave pair- 
ing interaction fixed to the crystal a-axis direction for all 
magnetization M _L a directions and the two parallel- 
spin states indicated by b = | tt) an( i c = | H) [8( . 

Subsequent measurements on a URhGe sample with 
RRR = 50 @ observed an anomalous high magnetic field 
H reentrant superconducting phase for iT||6 @, which 



was associated with a first-order metamagnetic transition 
and a topological Fermi surface Lifshitz transition 0, 0] . 
But the angular dependence of B C 2 within the ab plane in 
the low-field regime was found to be consistent with that 
of an ordinary superconductor with just effective mass 
anisotropy, at least within the relatively low resolution 
of the experiment designed for operation at much higher 
magnetic fields H [5]. At first sight, these results appear 
to be in contradiction with the earlier measurements of 
B c2 in URhGe Q. 

In this letter, we address this issue by calculating the 
full (8, (f>) dependencies of B c2 for the candidate p-wave 
polar/CBS state. We calculate the 8 dependencies within 
the a&-plane for the RRR = 21 and 50 crystals, and 
predict that under some conditions, URhGe could ex- 
hibit a non-monotonic _B C 2 (#,</>) curve with a peak at 
0° < 8* < 90° at fixed </>, providing a definitive bulk test 
of the orbital symmetry of the order parameter. Our 
method is applicable to all superconductors. 

Since T c in ferromagnetic superconductors such as 
URhGe is relatively low, we assume weak coupling for 
a clean homogeneous type-II parallel-spin p-wave super- 



conductor with effective Hamiltonian 12 



H = 
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where p — hk, k a = k ■ a, the orbital symmetry of 
the equal-spin pairing interaction is that of a p wave 
locked onto the orthorhombic crystal lattice & = e 3 
direction, and the spins indexed by a are quantized 
along the direction of B = (i H + M = V x A = 
B (sin 8 cos 4>, sin 8 sin </>, cos 8) . This model includes the 
spontaneous magnetization iWo of the ferromagnetic su- 
perconductor, which produced only a minor jump in B C 2 
near to T c in orthorhombic URhGe Here we assume 
an ellipsoidal Fermi surface with single-particle effective 



2 



masses mi along the orthogonal directions, respec- 
tively. These rtii are generally different from anisotropic 
Ginzburg-Landau (AGL) effective masses jljjj] . 

The Maxwell eq uation-preserving Klemm-Clem trans- 
it, 2(1 are exact at B C 2 in the AGL model, 



formations 



and were subsequently applied to a microscopic calcula- 
tion of B C 2 in (i-wave superconductors with mi = to 2 < 
Here the transformations allow us to calcu- 
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late the effects of a general ellipsoidal Fermi surface on 
B C 2 for a p-wave superconductor in the polar/CBS state, 
for which the order parameter anisotropy has a much 
stronger effect upon J3 c2 (0, <f>) than in those (i-wave cases. 

Here the anisotropic-scale transformation changes the 
Fermi surface shape from ellipsoidal to spherical 2(J 21 1 , 
allowing us to employ the Feynman technique used to 
calculate -B c2 (T) for an isotropic s-wave superconductor 
neglecting Pauli limiting [22j. We set pi — p^y/Wl, 
Ai = A'i^/Wl, Xi — x'i/^JWi, and Bi = B[j \JW~i, 

1 /3 

where m; = rrii/ (m\m2mz) . This transforma- 
tion changes B to B' (sin0'cos0', sin0'sin</>', cos9'), 
where sin0' = (m^^sin^/m^^), sin0' = 
[mi2(4')] 1 ^ 2 sm9/a(9, (/>), rai 2 (</>) = micos 2 ^ + m 2 sin 2 0, 



a(9,(j)) = [mi 2 (0) sin 2 9 + m 3 cos^0J ~' etc 

Next, we rotate the scaled coordinates so that the ro- 
tated induction B" points along the original pairing e 3 ' 
direction. Finally, we make an isotropic scale transforma- 
tion by l/a(9, 4>) [19j, |20|. The transformations have two 
overall effects: First, B -> Ba(9, (/>), modifying the slope 



,1/2 
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of b C 2 at T c , even for an s-wave superconductor Til- 23 1. 
Second, the rotation modifies V aa (k, k') to 



Vaa{k,k) = 3V o (fe3Cos0'-fc 2 sin0') 
x (fe 3 cos0' — fc 2 sin0'), 



(3) 



which differently alters b C 2(0,4>,t) from that of its slope 
at T c . We insert this rotated interaction into the lin- 
earized Gor'kov gap equation [2(j, project out the or- 
der parameter component, expand it in terms of the vor- 
tex harmonic oscillator states 12, 13j], and obtain a gen- 



eral recursion relation for the expansion coefficients a n , 
T n a n = 5sin 2 0'(/3„a„ +2 + f3„ 



r„ 

v (P.a) 



= -hit 



cos 2 0'a(f) 4 
d9 k sm9 k 



2^-2) > where 

(a) 



sin 



3 cos 2 
| sin 2 ! 



die 



ttTV ^ d9 k hm 3 9 k [ 
Jo 1 Jo 



xe 



p=0 



d£ e -«l«»l 

+ l)(n + 2)n! 
p!(p + 2)!(n-p)! : 



(4) 



(5) 



(6) 



t 



T/T c , T c = (27wo/7r)exp(-l/JV(0)Vb), JV(0) is 
the density of states per spin at the Fermi energy, u>q 



is a characteristic pairing cutoff frequency, 7 « 1.781 
is the exponential of Euler's constant, L n (z) are the 
Laguerre polynomials, and r\ = eB£ t 2 VpSm 2 9 k a(9, <fi). 
Equations (4-6) differ from those in Refs. [12,13] only 
in the general 0' in Eq. (4) and in the factor a(0, 4>) in rj. 
From the recursion relation above Eq. (4), S c2 (0, <fi, t) is 
implicitly obtained from the continued-fraction equation, 
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Usually, 4 or 5 iterations yield sufficient accuracy. 

In Fig. 1 (a), we assumed a spherical Fermi surface and 
plotted the reduced induction b c2 — 2eB c2 Vp / \2irT c ) 2 
versus t for 9 values increasing from 0° to 90° from top to 
bottom in increments of 10°. fr c2 (0°, t) and 6 c2 (90°, t) are 
those of the polar and CBS states [13J. 6 c2 (0, t) decreases 
monotonically with increasing 0, but is less sensitive to 
for 0^0° and especially for ~ 90° than in the case of 
ordinary effective mass anisotropy. As increases from 
0° to 90°, ~db C 2(9,t)/dt\t=i decreases monotonically by 
an overall factor of Since this slope variation is 

indistinguishable from that which could arise from ef- 
fective mass anisotropy, the same curves are rescaled by 
— db C 2/dt\t=i in Fig. 1(b). Order parameter anisotropy 
effects are easiest to identify for tel. 13l. Il4j. 




FIG. 1: (a) Plots of 6 c2 (0,i) = 2eB c2 v 2 F / \2nT c ) 2 for a spheri- 
cal Fermi surface with 9 increasing from 0° (top, polar state) 
to 90° (bottom, CBS state) in increments of 10°. (b) The 
same curves in Fig. 1(a) curves normalized by —db C 2/dt\t=i. 

At fixed t, b C 2(9,<p,t) for a polar/CBS p-wave su- 
perconductor with a general ellipsoidal Fermi surface 
only depends upon a(9, </>) and sin 2 0', 6 c2 (7r — 0, <fi, t) = 
b c2 {9, 4>, t), the entire 4> dependence of 6 c2 can be incorpo- 
rated into the Fermi surface anisotropy function [lj| ^ 



7 2 (^) = m 3 /m 12 (4>), 



(8) 
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and -db c2 (9,(t),t)/dt\ t=1 oc [3 sin 2 (9/ 7 2 (» +cos 2 6»]- 1 /2 
;13], suggesting 7 2 (</>) = 3 signals a crossover from order 
parameter to Fermi surface anisotropy asf^ 1 — . 

In Fig. 2, we plotted b C 2{6,t)/b C 2(0,t) for a variety of 
fixed 7 2 (<^>) values at t = 0, \. At lower t and as 7 2 (0) 
increases from 0.1 to 3, there is an increasing difference 
between b C 2(0,t) and the effective anisotropic mass form, 

Vg{0) = [cos 2 0/6 2 2 (O°) + sin 2 0/6 2 2 (9O°)]" 1/2 (9) 

fitted at each t. For fixed 7 2 (0) increasing above 3, an 
anomalous peak at 0° < 9* < 90° is clearly evident at t = 
and in the 3.7 and 5.9 curves at t = |. For 7 2 (</>) = 10, 
b c2 (9, \) only has a conventional maximum at 9 = 90°. 
The anomalous non-monotonicity is due to competing 
order parameter and Fermi surface anisotropy effects. 
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FIG. 2: (color online) Calculated b c2 {9,t)/b c2 (0,t) (solid) and 
fitted &®2 (#i t)/b c2 (0, 0), Eq. (9), (dashed) curves at constant 
7 2 (</>) values. The arrows indicate peak maxima at 9* points, 
(a) t = (b) t — 1/2. The inset is an enlargement of the 
80° < 9 < 90° region of the -y 2 (<f>) = 5.9 curves, with the 
indicated vertical scale points 1.2545 and 1.2549. 

We calculated the Fermi surface effective masses from 
the RRR = 21 URhGe crystal data [4]. In Fig. 3(a) we 
present the calculated b C 2(9, t)/b C 2(0, 0) in the ab plane 
(with B _L c) for different t values as functions of 9. 
The dashed lines represent fits to the corresponding fit- 
ted curve using Eq. (9). For low t, the order parame- 
ter anisotropy plays an integral role in the anisotropy of 
b C 2, while for t *~ 1, order parameter anisotropy is very 
small. Since the Fermi surface effective mass anisotropy 
is weaker in the ab plane than in the ac plane, our results 
differ substantially in this plane from those of Eq. (9). As 
noted above, in the be plane (9 — it/2), b c2 ((t)) oc 7(0). 
This is consistent with their data, since their b and c-axis 
data were both fit to the b c2 (t) of the CBS state 

We also performed similar calculations to make predic- 
tions for b C 2{9, (f>, t) in the RRR = 50 URhGe sample 5]. 
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FIG. 3: (color online) Calculated b c2 (6, t)/b c2 (0, 0) (solid) and 
fitted K2(8,t)/b c2 {Q, 0), Eq. (9), (dashed) curves for B _L c 
at various t values for the Fermi surface effective mass values 
obtained from experiment, (a) URhGe sample with RRR = 
21 Q. (b) URhGe sample with RRR = 50 @]. 
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FIG. 4: (color online) Logarithmic plot of j 2 (<f>) as a function 
of 9*, the peak angle in b C 2(0,t), at the indicated t values. 
Inset (a): Plot of the < 9* < 90° region versus log io [7 2 (0)] 
and t. Inset (b): Plot of b c2 (9, 0) /b c2 (0, 0) versus 9 near to 9* 
for 7 2 (</>) = 10 4 . The vertical scale runs from 46.5 to 47. 

In Fig. 3(b), the calculated b c2 (9, t)/b C 2{0, 0) and corre- 
spondingly fitted curves arc plotted in the ab plane as a 
function of 9 for various t, including t = 0.16, the low- 
est measurement value [5j|. From subsequent Shubnikov- 
de Haas measurements [6j, a strong B\\b appears to in- 
crease the pairing interaction strength Vo and decrease 
the relevant effective Fermi velocity [i[ of the heavy- 
electron ellipsoidal Fermi surface responsible for the pair- 
ing jol— [5Tlj| . By rotating a strong B a limited 4> range 
in the be plane, they concluded that the cyclotron mass 
to* = -y/mi2(90° + 0)to3 of the small spin-aligned heavy- 
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electron Fermi surface was independent of <f>. Although 
they could not determine j 2 without rotating B out- 
side the be plane, they provided strong justification for 
an ellipsoidal Fermi surface, and indicated that the Fermi 
wave vector kF decreases dramatically with increasing B, 
vanishing at 15.5 T, at which a topological Lifshitz tran- 
sition occurs [f|. We note that b C 2 and h c % differ greatly 
for these field strengths due to the large Afo||c pj. More 
important, if the order parameter in the reentrant phase 
maintains the polar/CBS form dramatic further in- 
creases in Vq and potentially in 7 2 (</>) would be expected 
as the metamagnetic transition is approached Q , and the 
angle between B and H would decrease dramatically 0] , 
yielding an observable anomalous peak in b C 2{9,t) as de- 
picted in Fig. 2. We propose that further measurements 
of h C 2(9 7 t -C 1) be made in the low-field and especially 
in the reentrant phase of URhGe, in which M changes 
rapidly from c to b, as this could provide a definitive 
identification of the symmetry of the order parameter 
and further support to a ferromagnetic mechanism for 
superconductivity in both phases of URhGe. 

In Fig. 4, we plotted log 10 [7 2 (</>)] versus 9* , the anoma- 
lous peak angle in b C 2(9, 0). Anomalous peaks appear for 
7 2 {4>) > 3 and disappear for 7 2 (</>) < X(t), where X(t) 
increases very rapidly with decreasing t for t < 0.15, as 
shown in inset (a). Inset (b) details the anomalous peak 
in b c2 (9,0) for 7 2 (» = 10 4 . 

Conventional peaks in 6®f(6*) occur only at either 
9 = 0° or 9 = 90°, but anomalous peaks only oc- 
cur for < 9* < 90°. However, since b c2 (9, tf>, t) = 
fo C 2(180° — 9, tfi, t), a second anomalous peak at 180° — 9* 
is reflection-symmetric in shape about 90° to that of the 
first one. When 9* is close to 90°, the magnitude of 
each anomalous peak is very small, but accurate mea- 
surements of this double peak could provide a definitive 
bulk test of the orbital symmetry of the order parameter. 

Ellipsoidal Fermi surface anisotropy can be included in 
accurate calculations of B c2 (9,(j),t) for superconductors 
of any orbital order parameter symmetry. If S^RuO^ 
were a chiral p x ± ip y state as purported (24[, B r ^{9, t) 
would be that of an SK state fixed to the x y pl ane 121. 
although B c2 (90°,0) appears Pauli-limited [la IM E3 ■ 
An approximation of our procedure was employed to fit 
the extremely Pauli-limited in-plane £> C 2(90°, <f>, t <C 1) of 
CeCu2Si2 [27|. An accurate calculation of B c2 (9,cj),t) at 
intermediate 9 values could be more definitive. Detailed 
B C 2(9, <t>, t) for the proposed /-wave forms for the C phase 
of UPt 3 could provide supporting information [28| . 

From analytic expressions for parallel-spin, p-wave su- 
perconductors with completely broken symmetry, we cal- 
culated B c2 (9,(j),t) with general ellipsoidal Fermi sur- 
face anisotropy. For Fermi surface anisotropy function 
7 2 (^>) > 3, the competing effects of order parameter and 
Fermi surface anisotropy lead to a novel, non-monotonic 
9 dependence that can provide a definitive test of order 
parameter symmetry in URhGe. Our technique can gen- 



erally provide accurate bulk tests of order parameter or- 
bital symmetry for unconventional superconductors, pro- 
vided that Vcr.cr' (k, k) and when necessary, Pauli-limiting 
interactions, are appropriately transformed. 
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